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Abstract 

I show that the strong coupling solution of the Kazakov-Migdal model with a 
general interaction potential V(3>) in D dimensions coincides at large N with that 
of the hermitean one-matrix model with the potential V(&): 



(2D - 1)V' = (D — 1)V' + D^(V') 2 + 4(1 - 2L>)$ 2 , 

whose solution is known. The proof is given for an even potential V(#) = V{— $) 
by solving loop equations. 
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1 Introduction 



Solvable matrix models are usually associated with D < 1 dimensional theories. Kazakov 



where the scalar field <& x is in the adjoint representation of the gauge group SU(N C ) and 
the link variable U^x) belongs to the gauge group, is solvable in the limit of large number 
of colors, N c , for D > 1. Since stringy phenomena are associated with the strong coupling 
phase (see for a review), one is interested in the solution of KMM at strong coupling. 

The original idea to solve KMM is based on the fact that the scalar field can be 
diagonalized by a (local) gauge transformation so that only 0(N C ) degrees of freedom 
are left and the saddle-point method is applicable as N c — > oo. Migdal proposed to 
solve the saddle point equation by the Riemann-Hilbert method and derived a master 
field equation to determine the N c — oo solution in the strong coupling phase. An 
explicit solution of this equation for the quadratic potential was found by Gross ||. A 
surprising property of the master field equation (not yet completely understood) is that 
it admits 0, || self-consistent scaling solutions with non-trivial critical indices. 

There are two other approaches to solving KMM at strong coupling. The first one 
is based on loop equations which were derived and solved for KMM with the quadratic 
potential by the author |J recovering the solution of Ref. [|]]. The second one, proposed 
by Boulatov 0, relies on the relation of KMM to a matrix model on the Bethe lattice, 
whose "naive" continuum limit is equivalent for D > 1 to a one-matrix model with the 
upside-down potential. 

In the present paper I extend the approach based on loop equations to the case of 
KMM with an arbitrary potential. On the one hand, this provides a way to solve the 
model at strong coupling which is an alternative to the Riemann-Hilbert method. On 
the other hand, this method is along the line of modern studies of matrix models of 2D 
quantum gravity by means of loop equations (for a review, see Ref. ||). 

The main results of this paper are as follows. The loop equations are drastically 
simplified in the strong coupling phase at N c = oo due to the fact that the averages 
of closed Wilson loops vanish except for the loops with vanishing minimal area. The 
resulting equations are satisfied for an even potential V(<&) = V(— <&) by the ansatz which 
reduces KMM in D dimensions to a hermitean one-matrix model with the potential V(<&): 



and Migdal [|T| have recently proposed that the model (KMM) defined by the partition 
function 




UdU,(x)U^e^<- V ^ + ^ 
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(V'($))2 + 4(1 - 2£>)$ 2 . 
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The solution of this model with an arbitrary potential is well-known. 
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In Section 2 I derive the exact loop equations for KMM with an arbitrary potential 
and any N c . In Section 3 I show how these equations are simplified in the strong coupling 
phase at N c = oo. In Section 4 I define the ansatz and obtain the exact solution of the 
N c — oo loop equations in the case of the even potential. Section 5 is devoted to the 
explanation of the form of the solution from the viewpoint of the large mass expansion. 
In Section 6 I discuss some properties of the solution and show how it can be compared 
to the one obtained by the Riemann-Hilbert method. Appendix A contains the details 
of derivation of the loop equations. In Appendix B I analyze the large mass expansion of 
the one-link correlator of the gauge fields at N c = oo. 



2 Loop equations for arbitrary potential 

The loop equations of KMM relate the closed adjoint Wilson loops 

W A (C) = ^(\trU(C)\ 2 -l)^ (2.1) 
to the open ones with the matter field attached at the ends: 

Gx(C xy ) = /i-tr^^C^^-C/t^))^ . (2.2) 

The loop equations result from the invariance of the measure in Eq. ( |1.1| ) under an 
arbitrary shift of $ and read 

/^^(y'i^mc^j^U^C^))) - £ G x (C {x+fl)x C xyJ 





U\C xy ))) (2.3) 



where the path C( x+fl ) x C xy on the l.h.s. is obtained by attaching the link (x, /z) to the path 
C xy at the end point x as is depicted in Fig. [p. The details of derivation are presented in 
Appendix A. I have omitted here and below additional contact terms which arise at finite 



N c due to the fact that <3> belongs to the adjoint representation, so that Eq. (pO|) is written 
for the hermitean matrices. This difference should disappear, however, as N c — > oo which 
can be easily proven for the even potential when (-^ tr V'(Q X )) vanishes. 

The path C xy on the r.h.s. of Eq. (|2.3|) is always closed due to the presence of the 
delta-function. The explicit equation for the vanishing contour C xx = at large N c , when 
the factorization holds, reads 

/ ^ V ^-2DG x( l) = El (2.4) 
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Fig. 1: The graphic representation for G\(C xy ) (a) and G \(C( x+ ^ x C xy ) (b) enter- 



ing Eq. ( |2.3| ). The empty circles represent Q x or <& x +fj, while the filled ones 
represent . The oriented solid lines represent the path-ordered products 
U{C xy ) and U{C^ x+ ^ x C xy ). The color indices are contracted according to the 
arrows. 



where 



A - $, 



X (Ga(0) + 1) 



with G A defined by Eq. ( [2.2|) . I have denoted the one-link average by 



(2.5) 



(2.6) 



since the r.h.s. does not depend on x and fi due to the invariance under translations by 
a multiple of the lattice spacing and/or rotations by a multiple of tt/2 on the lattice. 
The contour C\ encircles singularities of E w so that the integration over u on the l.h.s. of 
Eq. ( p.4[) plays the role of a projector picking up negative powers of A. 



3 Loop equations at large 7V C 

The loop equations for non-vanishing contours C xy ^ are drastically simplified at 
N c — oo in the strong coupling region where the closed adjoint Wilson loops (|2 . 1|) vanish 
except the contractable loops (i.e. those with vanishing minimal area A m i n (C) which are 
equivalent to C xx = due to the unitarity of i7's): 

W A (C) = 6 0Aminic) +o(±). (3.1) 
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While the averages of a new kind arise on the r.h.s. of Eq. ( |2.3| ) for C xx 7^ 0, they obey 
at N c = 00 the following analogue of Eq. ( |3.1| ) 



(U(C xx )-L-)±tr (u\C xx )-±-)\ = S , Amin(c) El + O (±) (3.2) 



\jV c V V XX '\-®jN c V V XXJ \ I U > A ™»^ A 

i.e. vanish for C xx 7^ 0. 

Hence, the strong coupling loop equation for C xy ^ at N c — 00 reads 

^i-tr (y\§ x )U{C xy ) J ^U^C xy ))^ - £ Gx(C (x+ , )x C xy ) = (3.3) 

independently of whether C xy is closed or open. 

Therefore, the r.h.s. of the loop equation in nonvanishing at N c = 00 only for C xy = 



(modulo backtrackings) when the proper equation is given by Eq. (|2.4j) . This property of 
the strong coupling loop equations at N c = 00 allows to find a simple solution. 



4 The strong coupling solution 



Let us solve the set (|2.4j) , ( p.3|) of the N c = 00 loop equations at strong coupling by the 
following ansatz in the case of the even potential 

^tr(F($ x )f/ M (x)$, + ^ (x) ^ = (itr(F(S <B )S a( A(S,))) (4.1) 

where F(<&) is arbitrary. The function A (a;) is analytic at ui = 0: 

00 

AM = ^A^. (4.2) 

k=0 

For G\(l) which is defined by Eq. (|2.6| ), Eq. ( [4.1| ) can be written as 

, to<gw K (4JJ) 



where the contour G\ encircles singularities of E w , i.e. the same as in Eq. (fZ 
The formula (|4.3|) extends to the general potential the one 



G\{1) = A Q \E\ (quadratic potential) (4.4) 
for the quadratic potential which is associated with || 

A (a;) = A (quadratic potential) , (4.5) 

i.e. Afc = for k > 1. 
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The substitution of the ansatz (|4.3| ) into Eq. ( |2.4| ) yields 

duj V'iuo 



E u = E{ (4.6) 



'ci 2ni A — uj 
where 

V'{uj) = V\u) - 2DujA(uj) (4.7) 

which coincides with the loop equation for the hermitean one-matrix model with the po- 
tential V. i 

The dependence of A (a;) on the potential V can be determined from Eq. (|3.3|) . The 
simplest way to do this is to consider Eq. ( |3.3|) in the case when C xy is just one link (x, /io) 
and to take the 1/ A 2 term of the 1/A expansion. The resulting equation reads explicitly 

" E (^tr(^_^(x-//)^ (x)$ s+w C/t o (x)^t( a; _ // ))) = o (U 
which reduces after the substitution of the ansatz (|4.3|) to 



I (n^ AM _ ! _ (2D _ 1)A 2 (a ^ ^ = o . ( 4. 9) 

JCi 2,711 \ 00 J 

This equation is satisfied provided that A(u) obeys the quadratic equation 

_ _L _ (2D - 1)A(«) = (4.10) 
w A(u;) 

for the generic potential 

oo 

vh = E^ 2 " ( 4 - n ) 

n=l 

and ^2 = ttiq. 

The solution to Eq. (|4.10|) reads 

AM = 1 2 (4.12) 

XM + ./«^ 2 + 4 (l-2D) 



which recovers the one || for the quadratic potential when 

V\uj) = 2m\ijj (quadratic potential) (4-13) 

and Eq. ( |4.5|) holds. The fact that the ansatz (|4.3|) with A[uS) given by Eq. Q4.12j ) satisfies 
Eq. ( |3.3|) means that it is indeed a solution providing Eq. (f4.6|) with 



2D-1 u 2D-l\ \ uj ) y ' V ; 



2 For a review, see Ref. 
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is satisfied. 

The one-cut solution to this equation for an arbitrary potential is well-known 



cluj V'(uj) 



(A-x)(A-y) 

'd 47ri A — uj \ {uj — x) {uj — y) 

where x and y are expressed via V by 



(4.15) 



duj 



V'{u>) 







do; 



o;y'(o;) 



Oi 2ttz ^1 {u - x ){u} - y) 



(4.16) 



One gets a; = — ?/ for the even potential. 

The formulas ( |4.3| ), ( [4.14| ), ([4.1 5|) and (|4.16|) completes the solution of the strong 
coupling loop equations of KMM at N c = oo. The saddle point value of $ s is totally 
determined to be ^-independent and is described, modulo a gauge transformation, by the 
spectral density 



P(A) 



2vr 2 



(x-\)(\- V ) 



for y < X < x 



t-X \| (x-t)(t-y) 
with support y < X < x. Properties of this solution are discussed in the next section. 



(4.17) 



5 Relation to the large mass expansion 



The peculiar form ( |4.3|) of the strong coupling solution at N c = oo can be understood 
in the framework of the large mass expansion. To this aim let us consider the one-link 
correlator of the gauge fields 



/ 1 , „M / 



(^tr ^wox. = j — — ^ ( ,::;\.\ '"' 1 (5-1 



J all e { 



where the averaging is only w.r.t. U while Q x and $ X + M play the role of external fields. t a 
(a = 1, . . . , AT? — 1) stand for generators of SU(N C ) which are normalized by 

tr t a t b = 5 ab . (5.2) 

As was proposed in Refs. [[3], fL0|| , the following formula holds at N c = oo: 

1 \ 00 1 

- tr {fUS^lA ) = £ A m - tr 1 ) . (5.3) 

It is shown in Appendix B how this formula can be obtained in the large mass expansion. 
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Eq. (|5.3|) allows to explain the solution of the previous section as follows. Let us 

by tr(t a F(§ x )) which gives, using the completeness 



multiply both sides of Eq. 
condition ( |A.4| ), 



i-tr (/•(<i. r )r<i>,.. / ,r ; 

^ tr (f(*,)*«A(*.)) + ^ tr (F(*,)) ^ tr - $,A($,)) (5.4) 

where A(<&) is defined by Eq. (|4.2|) . The second term on the r.h.s. which is due to the 
difference between the adjoint representation and the hermitean matrices vanishes for the 
even potential at A" c = oo. Hence, Eq. (|5.4j) for $ s and given by the saddle point 
matrix $s recovers Eq. (JO). On the other hand, the solution of the previous section 



allows to calculate A m in Eq. ( |5.3| ) as the coefficients on the expansion of ( |4.12| ) in 



6 Discussion 

The above strong coupling solution is realized at given D only in some region of the 
couplings t2fc's entering the potential ( EOT ). At D = it coincides with the well-known 



solution of the hermitean one-matrix model. At any D but t 2 ^ 0, t 2 k — for k > 1, the 
solution coincides with the one for the quadratic potential Q. For this reason I expect 
that it is realized in some region around this point. The condition is that Eq. (|4.16|) 



should yield real x and y and the spectral density (|4.17|) , which describes the distribution 
of eigenvalues of the saddle point matrix $5, should be positive. This is a restriction on 
the one-cut solution which is satisfied in some region of values of the couplings t^'s. It is 
well known that for the simple quartic potential 

V(u) =i 2 uj 2 + I 4 uj a (6.1) 

the one-cut spectral density is positive for ^1 < t\jY2. This is, however, a nontrivial 
restriction on the potential V since 

V\u) = D^{V'{u)f + — (D — l)V'{u) . (6.2) 

One more restriction on the one-cut solution is given by the requirement that the 
expession under the square root in Eq. (|4.14 ) must be positive for any u which belongs to 



the support of the spectral density. If this expression becomes negative for some values of 
t2fc's this simply means that the one-cut solution is not realized and one should look for 
a more sofisticated support (multi-cut solutions). This is quite standard for the large- A" 
phase transitions which occur at the values of couplings where the behavior of the spectral 
density changes. 

It is interesting to compare our solution with that obtained by the Riemann-Hilbert 
method. It is easy to identify the function T\(z) which determines the solution of Ref. 0] 
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with the one-link correlator 

r AW-l = /i?rtr(^— U— ^— E7+)\ (6.3) 

which is defined by the same average as in Eq. fl5.1| ). The proof is based on the following 
extension of the Migdal procedure. Let us define [[| 

Ca(^) = . 1 , i a *W) (6-4) 

where J stands for the Itzykson-Zuber integral 

= / ^e Wc "( WWt ) • (6.5) 
By a direct differentiation of Eq. (|6.5|) one gets 



Rewriting the l.h.s. via the spectral density, substituting for $> x and ^ x +^ the saddle point 
value $5 and remembering the definition of T\(z) 



T aW = 1+/^*», (6.7) 
J // — A 

one proves Eq. ( |6.3|) . 

From Eq. (|6.6|) it is easy to see alternatively that N c = oo 

where the average is w.r.t. the same measure as in Eq. ( |1 . 1|) . Therefore, we get asymp- 
totically 

T x {z) = 1 + — + + • • • as z oo (6.9) 

z z z 

so that our G\(l) is to be compared with the 0(z~ 2 ) term in the expansion of T\{z). 
It would be very interesting to calculate exactly the correlator on the r.h.s. of Eq. ( |6.8|) . 
which should be manifestly symmetric w.r.t. A and z, by solving the loop equations in 
order to compare with T\(z) of Ref. ||. 

It is worth mentioning that in finding the solution I did not use the fact that Q x can 
be diagonalized. For this reason a solution analogous to that of this paper exists for 
the adjoint fermion model |l ij which reduces to the complex one-matrix model [0]. This 
result, as well as an analysis of the D < 1 case, will be published elsewhere. The most 
interesting question is whether the solution of this paper admits a continuum limit for 
D > 1. 
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Appendix A Derivation of the loop equations 



Let us consider an equation which results from the invariance of the measure over $ in 
the open-loop average fl2.2| ) under an infinitesimal shift 

— + & (A.l) 

of $ x at the given site x with £ x being an infinitesimal hermitean matrix. For KMM 
one should impose tr^ = in order for the shifted matrix to belong to the adjoint 
representation of SU(N C ). Since this condition should be inessential as iV c — > oo, I derive 



loop equations for the hermitean model, defined by the partition function ( |1 . 1| ) with the 
integration going over arbitrary hermitean matrices $> x . £ x in Eq. ( A.l|) is then arbitrary 
hermitean. 

It is convenient to introduce N% generators 

[t A ] tl = fa, mj) (a.2) 

with t a (a = 1, . . . , N%—1) being the standard generators of SU(N C ). The generators (|A.2|) 
obey the following normalization 

-L t rt A t B = 5 AB (A.3) 

and completeness condition 

[t A ] ij [t A ]ki = Hc8ah j - (A.4) 
An arbitrary N c x N c hermitean matrix $ can be represented as 

$ = t A< t> A where $ A = (— tr $, — trt a $) (A.5) 

with $° = tr $ vanishing if $ is taken in the adjoint representation of SU(N C ). 

To derive the loop equation I apply a trick similar to that used in deriving loop 
equations of QCD ||. Let us consider the loop average 

(^±t^t A U(C xy )^^U^C xy ))^=0, (A.6) 

where the averaging is taken with the same measure as in Eq. fll.lQ , which vanishes due 
to the gauge invariance. Performing the shift ( |A.1|) of <& x , using the invariance of the 
measure and calculating d/d§ B (x), one gets 

(w tr i tBv '^) w tr ( tAu ( c ^j^ ut{Cxy] 

\ c c y i 

- £ (^tr(t B C/ M (a;)^C/t(x))^tr^C/(a,) ) ^[/t(^)^ 

= *** (jrz tT ( tAU{C ^X^/x^ y Ut{C ^)) ■ (A ' 7) 



10 



The l.h.s. of this equation results from the variation of the action while the r.h.s. represents 
the commutator term resulting from the variation of the integrand. 

The averaging over the gauge group picks up two nonvanising invariant equations for 
the hermitean matrices. The first one can be obtained contracting Eq. 
(A, B = 0, . . . , N% — 1) while the second one is given by the A, B = component. 

The first equation reads 

D I 1 1 

where the contour C[ x+ ^ x C xy is obtained by attaching the link (x, /i) to the path C xy at 
the end point x as is depicted in Fig. |T|. Using the definition ( j2.2| ), this equation can be 
written finally in the form ( |2.3| ) . 

The second equation which is given by the A, B = component of Eq. ( |A.7| ) reads 
/J_t r y'($J— tr— \- T /— tr$, +it — tr — \ 

In the large- iV c limit when the factorization holds, Eq. (|A.9| ) is automatically satisfied as 
a consequence of the C(A _1 ) term in Eq. ( |A.8| ). 



O) by 5 



AB 



Appendix B The one-link correlator at large iV c 

Eq. (|5.3| ) can be derived for $ given by the saddle point (master field) configuration 
analyzing the large mass expansion which allows to calculate the one-link correlator (|5.1|) 
in the strong coupling phase (i.e. before an expected large- iV c phase transition). To 
calculate it, let us expand the numerator in powers of $ as is depicted in Fig. ||. 

The idea is now not to calculate the complicated integral over U^ix) in Eq. (|5.1|) 
but rather substitute for $ x and the saddle point value $5 which is determined 

by the future integration over <3> according to Eq. ( |1 . 1| ) . Notice, that the integrals over 
$ x and Qx+ft are independent to each order of the large mass expansion. Therefore, one 
substitutes 

(K)s«)s = K—6 ab (B.l) 
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t a 



Fig. 2: The graphic representation of the large mass expansion of the one- link corre- 
lator (5.1). The right filled circles represent ^ x +^i while the left ones represent 



<&x- The empty circle represents t a . 



and 
where 



ffls(0 = (B.2) 
K=^tr$| (B.3) 



and ( |B.2[ ) vanishes due to the gauge invariance. 

After the use of Eqs. ( |B.1| ), ( |B.2| ) and the completeness condition ( |A.4| ) the color 
indices are contracted in such a way that all the matrices U^(x) disappear due to the 
unitarity. The simplest contractions are depicted in Fig. |3|. Only the connected diagrams 
should be taken into account due the presence of the denominator. The diagrams of the 
type of Fig. are the only ones which emerge for the quadratic potential. They always 
result in Ao- The diagrams of the type Fig. |b appear when the interaction is present. 
They result in A m with m > 1. 

Finally, let us notice that this structure of A m is not spoiled by the fact that Q x enters 
2D links emanating from the point x. This affect only combinatorics making K to be 
D- dependent. 
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Fig. 3: The contraction of color indices for diagrams of Fig. |2[ 

a) The diagram results in the contribution to Ao- 

b) The diagram results in the contribution to Ai. 
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